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Solid-state qubits hold the promise to achieve unmatched combination of sensitivity and spatial 
resolution. To achieve their potential, the qubits need however to be shielded from the deleterious 
effects of the environment. While dynamical decoupling techniques can improve the coherence 
time, they impose a compromise between sensitivity and bandwidth, since to higher decoupling 
power correspond higher frequencies of the field to be measured. Moreover, the performance of 
pulse sequences is ultimately limited by control bounds and errors. Here we analyze a versatile 
alternative based on continuous driving. We find that continuous dynamical decoupling schemes 
can be used for AC magnetometry, providing similar frequency constraints on the AC field and 
improved sensitivity for some noise regimes. In addition, the flexibility of phase and amplitude 
modulation could yield superior robustness to driving errors and a better adaptability to external 
experimental scenarios. 



Solid-state qubits have enierged as promising quan- 
tum sensors, as they can be fabricated in small volumes 
and brought close to the field to be detected. Notably, 
Nitrogen- Vacancy (NV) centers in nano-crystals of dia- 
mond [1] have been applied for high sensitivity detection 
of magnetic [2-4] and electric fields [5] and could be used 
either as nano-scale scanning tips [6] or even in- vivo due 
their small dimensions and low cytotoxicity [7]. Unfor- 
tunately, solid-state qubits are also sensitive probes of 
their environment [8, 9] and this leads to rapid signal de- 
cay, which limits the sensor interrogation time and thus 
its sensitivity. Dynamical decoupling (DD) methods [10- 
14] have been adopted to prolong the coherence time of 
the sensor qubits [2, 8, 15, 16]. Although DD techniques 
prevent measuring constant, DC fields, they provide su- 
perior sensitivity to oscillating AC fields, as they can in- 
crease the sensor coherence time by orders of magnitude. 
The sensitivity is maximized by carefully matching the 
decoupling period to the AC field; conversely, one can 
study the response of a decoupling scheme to fields of 
various frequencies, thus mapping out their bandwidth. 
Still, the refocusing power of pulsed DD techniques is ul- 
timately limited by pulse errors and bounds in the driv- 
ing power. Here we investigate an alternative strategy, 
based on continuous dynamical decoupling (CoDD), that 
has the potential to overcome these limitations. 

We consider the problem of measuring a small ex- 
ternal field, coupled to the sensor by a Hamiltonian: 
T-Lf, = jb{t)Sz, where Sz is the spin operator of the quan- 
tum sensor. For example, b{t) can be an external mag- 
netic field and 7 the spin's gyromagnetic ratio. The figure 
of merit for a quantum sensor is the smallest field Sbmin 
that can be read out during a total time t, that is, the 
sensitivity 77 = 6bmi7iVi- We use this metric to compare 
pulsed and continuous DD schemes and show how CoDD 
can offer an advantage for some noise regimes. 

The principle of DD schemes rests on the spin echo 
sequence, which refocuses unwanted phase accumulation 
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FIG. 1. Pulse sequences for four AC magnetometry schemes: 
PDD (P), constant driving (S), RE with optimal frequency 
(R^''') and spin-locking (S). Blue boxes represent microwave 
driving, with phase (x and y) as indicated. 

due to a slow bath by reversing the system evolution 
with control pulses. More complex DD sequences can 
in principle extend the coherence time indefinitely, by 
increasing the number of pulses. In practice, however, 
a large number of imperfect, finite-width pulses pro- 
vokes the accumulation of error and degrades DD per- 
formance [13, 17, 18]. CoDD has been first introduced 
in the context of NMR to mitigate pulse errors [19, 20] 
and it has then lead to many schemes, such as compos- 
ite pulses [21, 22], dynamically corrected gates [23] and 
optimized modulations [24] . In general, phase and ampli- 
tude modulation of the continuous driving allows great 
fiexibility and CoDD can achieve high decoupling power. 
Here we consider only two schemes, constant continuous 
driving (C) and Rotary Echo (RE) [25-27], as their pe- 
riodicity allows an easier use for AC magnetometry (see 
Fig. 1); we will compare these schemes to the simplest 
pulsed DD scheme, period dynamical decoupling (PDD). 

As an example, we compute the signal and sensitivity 
of AC magnetometry under RE, but similar derivations 
apply for the other schemes. The RE sequence consists of 
a continuous on-rcsonance driving field of constant am- 
plitude n and phase inverted at periodic intervals (see 
Fig. 1). RE is parametrized by the angle — ilT/2, 
where T is the sequence period. While RE is usually em- 
ployed to refocus errors in the driving field, for — 27r/c 
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the sequence also refocus dephasing noise, with perfor- 
mance depending on both k and the Rabi frequency. We 
consider the evolution of a sensor qubit under a sequence 
of 27rfc-RE and in the presence of an external AC mag- 
netic field of frequency uj whose magnitude b is to be 
sensed: 



(1) 



where E>W{t) is the square wave of period T = Airk/Q,. In 
the toggling frame of the driving field, the Hamiltonian 
becomes 
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We consider only the cases where (p = and ujT = 2rmT, 
with m an odd integer, since as we show below this yields 
good sensitivities. Under this assumption 'H(t) is peri- 
odic and for small fields b the evolution operator can be 
well approximated from a first order average Hamiltonian 
over the period T, "H « ^ 'H{t)dt = jbSy. 

If m = 1, we define oJiow = H, which, for a fixed 
fi, is easily adjustable by changing the echo angle 27rfc. 
Setting instead m = {2k — 1), we define ojopt = ^^^2k^^ : 
which yields b = 4&fc/[7r(4fc — 1)] and attains the best 
sensitivity of the method. The sensitivity, obtained as 

7]{t) ~ lim T^i-jVt, where S is the signal and AS its shot- 
fc— >-o I "aF I 

noise limited uncertainty, depends on b, that is, on the 
averaging of the AC field over the sequence period due 
to the DD modulation. We compare the performance 
of both 27rfc-RE schemes to FDD (optimum cj — 27r/i, 
if — 7r/2) and a constant modulation with uj ^ Q (see 
Fig. 1). We obtain for the schemes considered: 
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where ry 



= , with C a parameter capturing inef- 



ficiencies in the sensor readout [2]. Here Rk labels a 
2fc7r-RE scheme, P the FDD scheme and C the constant 
modulation (see Figure 1). A fourth operating scheme 
can be obtained by a "spin-locking" sequence [28] , where 
the spin is first rotated to the transverse plane before 
applying a driving field in the same direction; choosing 
(p = and UJ — fl yields the same sensitivity as for the 
constant modulation, rjs = rjc, even when the driving 
phase is inverted periodically. We note that if the phase 
if of the AC field is not optimized, the sensitivities are 
reduced by a factor $(</j), with $p = = csc{ip) and 
^Rk = sec ip. If in addition the phase of the AC field can- 
not be fixed, ^{p) = \/2 when considering the average 
signal over many realizations. 

These ideal sensitivities are degraded in the presence 
of noise and whenever the frequency of the AC field is 
not matched to the DD period. In the following we an- 
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FIG. 2. Bandwidth for AC magnetometry. We plot the weight 
functions WiiS) that scale the phase acquired during DD mag- 
netometry for AC fields of frequency lo. Left: we plot Wiio) 
and the envelop of its passband decay for FDD (blue dotted) , 
RE (fc = 1, red, thick) and constant driving (green, thin line) 
for n = 2 cycles, expressing the frequency in terms of the se- 
quence period. In the inset: we compare the main peak for 
n = \ (red, thick) and n = 10 (gray) for RE (k — 1) showing 
the reduction in bandwidth. Right: we compare W(aj) for 
continuous driving (green) and for RE with k = 1 (red, thick) 
and = 4 (gray, dashed), plotting as a function of uu in units 
of the Rabi frequency fi. 

alyze these two contributions, showing that they lead to 
a sensitivity t] ^ ri'D{t) /W{uj), where 2?(t) describes the 
decay under DD sequences and W{uj) is the reduction 
in the accumulated phase when the field frequency cj is 
suboptimal. 

Optimal sensitivities are obtained by carefully match- 
ing the period of the DD schemes to the oscillating field. 
In practice, however, when field frequencies are either un- 
known or known to a finite precision, it is of relevance to 
determine the bandwidth of the scheme and the deviation 
from optimum sensitivities. We estimate the bandwidth 
by calculating the phase accumulated by the sensor over 
the total interrogation time t = nT, Bt = /J b{t)f{t)dt, 
and examining the frequency dependence of its absolute 
value. For FDD, the filter function is fp{t) = SWp{t), 
the square wave with the period of the modulation. For 
continuous driving schemes such as RE and Rabi, f{t) 
is the strength of the toggling frame Hamiltonian. In 
particular, fn^ (t) = SW(i) sin(ilt) yielding the weight 
function Wp^iu) = \BpJu;)\/\Bn,iu;''n\-- 
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Wr^ has peaks {pass-bands) at a; = 27r(2(fe + p) — 1)/T, 
where p is an integer satisfying p > 1— fc. The lowest pass- 
band occurs for p = 1 — k, corresponding to wjou, = ^l/2k. 
The strongest peak is for p = at uiopt- Subsequent 
periodic peaks are attenuated from the symmetry point 
Lj = VL a.s ^ l^2_Q2| -The FWHM of the optimum peak in 

Wr,^{uj) decays as « where 7.58 « FWHM of the 
squared sine function, a result common to the other DD 
schemes. 

A similar calculation for the accumulated phase dur- 
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FIG. 3. Sensitivity for AC magnetometry. We compare the 
magnetic field sensitivity of a single NV center for FDD (left) 
and RE (fc = 1 center; k — 4 right). We assumed T2 — 500/is 
under OU noise (comparable to a ^'^C bath), yielding a decay 
cxe~"^ "^2)^ and a single readout with C = 0.03. A larger 
number of refocusing cycles (with shorter periods) achieves 
better sensitivity but can only detect higher frequencies, as 
shown by the color of the curves (right bar, MHz). 



can be calculated from the average of the superoperator 
(S(f)) = (Te~*-'o ^'^^ ), where we indicate by a hat the 



superoperators A = A(E) t — i- ^ A and T is the time or- 
dering operator. In turns, this can be approximated by 
the cumulants, « exp [— (ifi + K2 + ■ ■ ■)t], with 

the first cumulant Ki = and the second given by 

Jo Jo 
where the cumulant average is 

{H{h)H{t2))c = T{H{h)H{t2)) - {Hih)){H{t2)). 

In the toggling frame of the driving field, the 
stochastic Hamiltonian is Ti{t) = S{t)N{t) = 
6{t) [cos{ftt)(Tz + sm{flt)(Ty]. Then the second cumu- 
lant for n cycles is K2 = nA + OJ2^=i{i^~' [30], with 

4/gTr j 

Gj — and 



ing a PDD sequence indicates the existence of peaks at 
w = rmi/T, with m odd, whose intensity decays as 1/m. 
This slower decay than for the RE pass-bands could be 
beneficial if the goal is to detect fields of imknown fre- 
quencies. On the other hand, AC magnetometry un- 
der continuous driving or spin locking could be used 
for frequency-selective detection because M^c(w) has a 
unique peak at w = 2'k/T with FWHM on the same or- 
der of that for RE. A comparison of the different weight 
functions is depicted in Fig. 2. We note that while W{lo) 
describes the poor performance of DD schemes at detect- 
ing AC fields with unmatched frequencies, this property 
could in turn be used for frequency-selective measure- 
ments and even spectroscopy, by scanning the sequence 
period. While constant driving provides the best selec- 
tivity (canceling out higher octaves), RE provides more 
flexibility by changing both the period time and the angle 
27rfc, which allows more uniform noise cancellation. 

The refocusing power of RE can surpass that of pulsed 
schemes. Consider for example a noise with long cor- 
relation time Tc'. In this limit, the signals decays as 

-(T,at?/n\ ^i^i^ ^ _3^^ US- 
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ing a similar derivation [29, 30] 
PDD sequence is instead {Sp{t)) 
T\p = The sensitivities in Eq. (3) are further lim- 
ited by the signal decay T>{t) under the DD sequences. 
The achievable sensitivity is then a compromise between 
the refocusing power of the sequence used and the fre- 
quency that it allows detecting (Fig. 3). While the de- 
cay for pulsed DD has been widely studied, evolution 
under continuous DD is more complex [31]. We can 
estimate the RE decay to first leading order using a 
cumulant expansion [29, 30]. We assume a stochastic 
Hamiltonian, nit) = nSWkit)ax + S{t)a:,, where S{t) 
is an Ornstein-Uhlenbek noise with zero mean and au- 
tocorrelation function G{t) = a'^e~^, with a the dis- 
persion and Tc the correlation time. The signal decay 
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The cumulant can be written as 

a + P 



K2 = "-^Sl- 
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(see appendix for explicit expressions), yielding the signal 



1 + e""(cosh(7) + - sinh(7)) 
7 



Numerical simulations match well with these approxi- 
mate analytical results. 

The longer coherence time under the RE sequence can 
be exploited either to reach a better sensitivity for a 
given frequency or to measure lower frequency fields at a 
given sensitivity, as shown in figure 4. The achievable im- 
provement depends on the effective coherence time ratio, 
r = T2R/T2E, obtained from the two schemes. Because of 
the improved refocusing of RE with respect to PDD, the 
sensitivity can be improved for some noise regimes. In 
addition, RE-AC magnetometry provides the fiexibility 
of using larger angles (larger k) to allow for longer in- 
terrogation times (Figure 3) at lower frequencies, which 
could be beneficial in practical cases in combination with 
repeated readout schemes [26, 32]. 

We remark that besides the decay functions obtained 
above in the presence of dephasing noise, other sources 
of decay can arise from imperfect pulses or fiuctuations 
in the driving power. To this effect, RE provides a 
good protection against slow fiuctuation in the driving 
power [25, 26] and it is thus expected to achieve much 
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FIG. 4. Sensitivity for AC magnetometry. We compare 
the achievable sensitivity for constant driving (green, dash- 
dotted), FDD of n = 50 echoes (blue, dotted) and RE (27r- 
RE, red, achieving the same sensitivity of FDD at a lower 
frequency and Stt-RE, black, achieving better sensitivity than 
FDD at the same frequency). We assumed T2 = 500^s under 
OU noise, yielding a super-exponential decay oc "^2)^ 
and a single readout with C = 0.03. The decay of the con- 
stant (Rabi) driving was calculated following Ref. [31] for long 
Tc- The dashed, thin lines correspond to the ideal limit with 
no driving or pulse errors. 



better overall sensitivities than a continuous driving. In 
conclusion, we analyzed a novel scheme for AC magne- 
tometry based on continuous dynamical decoupHng and 
compared its performance to pulsed DD schemes. While 
we focused on the simplest DD sequences, we note that 
more complex driving, such as composite pulses [22, 26], 
could achieve even better refocusing of driving field in- 
stability and inhomogeneity while still providing compa- 
rable sensitivity. We further analyzed the response of AC 
magnetometry to fields of unknown frequencies, finding 
that some CoDD schemes (such as continuous driving or 
spin locking with alternating phases) are advantageous 
for spectroscopy. The sensitivity is ultimately limited not 
only by the theoretically achievable coherence time, but 
also by pulse errors or fluctuations in the driving field. 
While a full comparison of the limits due to imperfection 
in the control fields is beyond the scope of this work, the 
flexibility of CoDD schemes in modulating both phase 
and amplitude of the driving fleld can provide practical 
advantages, yielding a better compromise between the 
DD refocusing power and the frequencies of the fleld to 
be measured. 
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Appendix A: Cumulant 

We can calculate the time (ensemble) average of a time-ordered exponential operator by means of a cumulant 
expansion The first cumulant is zero since we assume a zero-average as zero. The second cumulant for the RE 
sequence is given by Eq. 5 with 
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